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Abstract 

The current normalization of the cross section of inclusive high-pr particle produc- 
tion in deuteron-gold collisions measured at RHIC relies on Glauber model calculations 
for the inelastic dAu cross section. These calculations should be corrected for diffrac- 
tion. Moreover, they miss the Gribov's inelastic shadowing which makes nuclei more 
transparent (color transparency) and reduce the inelastic cross section. The magni- 
tude of this effect rises with energy and one may anticipate it to affect dramatically the 
normalization of the RHIC data. We evaluate the inelastic shadowing corrections em- 
ploying the light-cone dipole formalism which effectively sums up multiple interactions 
in all orders. We found a rather modest correction factor for the current normalization 
of data for high-pr hadron production in d — Au collisions. The results of experiments 
insensitive to diffraction (PHENIX, PHOBOS) should be renormalized by about 20% 
down, while those which include diffraction (STAR), by only 10%. In spite of small- 
ness of the correction it eliminates the Cronin enhancement in the PHENIX data for 
pions. The largest theoretical uncertainty comes from the part of inelastic shadowing 
which is related to diffractive gluon radiation, or gluon shadowing. Our estimate is 
adjusted to data for the triple-Pomeron coupling and is small, however, other models 
do not have such a restrictions and predict much stronger gluon shadowing. Thus, 
one arrives at quite diverse predictions for the correction factor which may be even 
as small as K = 0.65. Therefore, one should admit that the current data for high-py 
hadron production in dAu collisions at RHIC cannot exclude in a model independent 
way a possibility of initial state suppression proposed by Kharzeev-Levin-McLerran. 
To settle this uncertainty one should directly measure of the inelastic d — Au cross 
sections at RHIC. Also collisions with a tagged spectator nucleon may serve as a sensi- 
tive probe for nuclear transparency and inelastic shadowing. We found an illuminating 
quantum-mechanical effect: the nucleus acts like a lens focusing spectators into a very 
narrow cone. 

* Based on lectures given by the author at Workshop on High pt Correlations at RHIC, Columbia 
University, May-June, 2003. 
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1 Introduction 



Recent data for high-p T hadron production in deuteron-gold collisions at y/s = 200 GeV 
at RHIC [1-3] demonstrate importance of these measurements for proper interpretation of 
data from heavy ion collisions. The observed nuclear effects at high-pr are pretty weak, 
the enhancement (Cronin effect) measured for pions by PHENIX is only about 10 — 20%, 
in accordance with expectation of [4] and with somewhat larger effect found in [5], while a 
suppression, rather than enhancement was predicted in [6]. To discriminate between these 
predictions the data should have at least few percent accuracy. 

In this notes we draw attention to the fact that only the shape of p^-distribution was mea- 
sured experimentally, while the normalization of the data is based on theoretical calculations 
which are not correct. Therefore, the reported results of deuteron-gold measurements [1-3] 
may be altered by more appropriate calculations. 

The nucleus to nucleon ratio demonstrating the well known Cronin effect [7] is defined 

as, 

da hA /d 2 p T 

Ra/n{Pt) = AdoM/dtpr • (1) 

At large pt, of the order of few GeV this ratio exceeds one, but eventually approaches one 
at very high p T as is expected according to /c T -factorization (it may even drop below one 
due to the EMC effect at large Bjorken x). 

Absolute values of the high-p^ nuclear cross sections are difficult to measure at RHIC, 
only the fraction of the total inelastic cross section, dN hA / d 2 pT is known. Then, one has to 
normalize it multiplying the fraction by the total inelastic cross section, 

o hA dN hA /d 2 p T 1 dN hA /d 2 p T 

n A /N{PT) A(J NN dNh N/ d 2 pT dN hN/ d 2 pT > W 



where 

N eM = A^ x , (3) 



in 



In some experiments the denominator in (1), da hN / d 2 px was directly measured or borrowed 
from other measurements, otherwise it should be corrected for diffractive dissociation of the 
colliding protons which possesses a large rapidity gap and escapes detection. In what follows 
we assume that the denominator in (1), da hN /d 2 pT was directly measured (see, however, 
discussion in Sect. 1.2.2) and concentrate on nuclear effects, i.e. the inelastic nuclear cross 
section a^ n A which was calculated in [1-3] in an oversimplified approach. 



1.1 Number of collisions: who is actually colliding? 

The Glauber approach is a model for the elastic hadron-nucleus amplitude. It is demon- 
strated in Appendix A how to calculate inelastic and quasielastic cross section using unitarity 
and completeness. The model does not say anything about exclusive channels of inelastic 
interaction. One can formally expand the Glauber exponential, and it looks like a series 
corresponding to different numbers of inelastic collisions of the same hadron and with the 
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same inelastic cross section. However, a high-energy hadron cannot interact inelastically 
many times, since the very first inelastic collision breaks down coherence between the con- 
stituents of the hadron. It takes a long time proportional to the energy to produce a leading 
hadron in final state. 

The cross section of inelastic hadron- nucleus collision, a^ 4 , is related to the probability 
for the incoming hadron to get the very first inelastic collision, usually on the nuclear 
surface. This is why o\, A oc A 2 / 3 . Since the process is fully inclusive, subsequent final state 
interactions do not affect the cross section due to completeness. 

N co u defined via expansion of the Glauber exponential term should not be treated as 
multiple sequential interactions of the projectile hadron (like expansion of the exponential 
describing the time dependence of particle decay does not mean that the particle can decay 
many times). After the first inelastic interaction the debris of the projectile hadron keep 
traveling through the nucleus, but their interactions apparently have little to do with the 
properties of the incoming hadron and its inelastic cross section. Formally, one can relate 
Ncoii to the mean number of the Pomerons which undergo unitarity cuts. The Abramovsky- 
Gribov-Kancheli (AGK) cutting rules [8] which are not proven in QCD, assume that these 
cuts have the same eikonal weights as given by the Glauber model. In this approach multiple 
interactions are not sequential (planar), but occur in parallel, i.e. they allow a simultaneous 
unitarity cut. In terms of the light-cone approach multiple interactions correspond to higher 
Fock states in the projectile hadron. The constituents of this states propagate through the 
nucleus and experience their first inelastic interactions independently of each other. The 
probability of such multiple interactions has little to do with the properties of the low Fock 
states which dominate the hadron-nucleon cross section, therefore it should not be expressed 
as a power of . At any rate, whether the AGK weights are correct or not, it is clear 
N co u cannot be treated as sequential interactions of the projectile hadron. 

The first inelastic collision of the incoming hadron is a soft color-exchange interaction. 
The projectile partons do not alter either their number (for a given Fock state), or their 
longitudinal momenta, but the whole system of partons acquire a color. Therefore, the 
remnants of the hadron turn out to be color connected to the remnants of the target. Then 
new partons are produced from vacuum (e.g. via the Schwinger mechanism) aiming to neu- 
tralize the color of the projectile partons. Their momenta are much smaller than of the 
projectile partons. Such an excited and color neutral partonic system keeps propagating 
through the medium and experiencing new soft color exchange interactions similar to ordi- 
nary hadrons. The corresponding cross section is subject to color screening and is controlled 
by the transverse, rather than longitudinal size of the system. 

From the practical point of view, there is nothing wrong in using N co u as a multiplication 
factor for hard reactions, since within the Glauber model it is proportional to the nuclear 
thickness function T^(b), i.e. to the number of opportunities for a parton to perform 
a hard process. Indeed, the projectile high energy partons participate in hard reactions 
independently of the accompanying partons, since color screening plays no role for a hard 
interaction. Moreover, naively, one may expect that this factor T% (b) (Tab i n the case of AB 
collision) is all one needs to normalize a hard process, and this normalization is independent 
of the soft cross section cr^ . However, N co u is defined for events where inelastic collision 
did happen. Therefore, it must be properly normalized by the probability for the incoming 
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hadron to make inelastic interaction at the given impact parameter, 

N (h\ - in A \ ' (A\ 

CM[) W 

Averaging this expression over inelastic collisions at different impact parameters one indeed 
arrives to the expression Eq. (3). 



1.2 Correcting data for Rd-Au 

The current analyses of RHIC data [1-3] calculate N coU in Glauber Monte Carlo model 
assuming af n N = 41 — 42 mb. In these notes we challenge these calculations and show that 
the published results for d — Au collisions are subject to important corrections and the 
conclusions are model dependent. 

There are two major corrections to be done to the Cronin ratio Eq. (2) measured at 
RHIC. We combine them in a correction factor K, 

RdA (Pt) = RT°{Pt) xK, (5) 

where 

K = K Gr x Kqi . (6) 

Here K Gr is the correction related to Gribov's inelastic shadowing missed in Glauber model 
calculations. It is introduced in Sect. 1.2.1 and calculated throughout the paper. 

Even within the Glauber model the calculations performed in [1-3] should be corrected 
by a factor Kq\. It originates from a more accurate treatment of the inelastic NN cross 
section which should correspond to the class of events selected for the analysis, as is explained 
in Sect. 1.2.2. This correction is calculated in Sect. 2. 

There is an additional correction which should be included into (6) if one needs to 
compare with theoretical predictions for the Cronin effect for pA collisions. It is related to 
the fact that the deuteron is a nucleus and is also subject to the Cronin effect. Therefore 
high-p-r enhancement in d — A must be somewhat stronger than in p — A collisions. This 
correction is evaluated in Sect. 7 and found rather small. 



1.2.1 Inelastic shadowing. 

It is known that Gribov's inelastic corrections [9] to the Glauber approximation make nuclear 
matter more transparent and reduce the hadron-nucleus cross sections compared to the 
Glauber model. This effect steeply rises with energy, as one can see from the example 
depicted in Fig. 1 for the total neutron-lead cross section measured and calculated in [10]. 
Apparently, the Glauber model overestimates the cross section, and the deviation rises with 
energy. Without a good theoretical input one cannot predict what will happen at the energy 
of RHIC which is 100 times higher than in fixed-target experiments at Fermilab. This is a 
serious challenge for the theory to calculate the inelastic dA cross section at these energies, 
and the results apparently will be model dependent. However, it is certain that the sign of 
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Figure 1: Data and calculations [10] for the total neutron-lead cross sec- 
tion as function of energy. The dashed curve corresponds to the Glauber 
model, while the solid curve is corrected for Gribov's inelastic shadowing. 



the correction remains negative and it can only rise with energy, i.e. cannot be smaller than 
what is shown in Fig. 1 for low energies. 

Our own estimates summarized in Table 1 give a moderate reduction, about 20%. The 
weakness of the effect is based on a proper treatments of diffraction and is fixed by data on 
large mass diffractive dissociation of protons [11]. At the same time, many models predict 
quite a strong gluon shadowing even at high virtualities. Naturally, this effect should not be 
weaker in soft NN interactions. Then it may lead to a stronger suppression of the inelastic 
dA cross section than we found, as it is discussed in Sect. 6.2. 

Note that although inelastic shadowing makes nuclear medium more transparent, the 
mean number of collisions increases according to (3). It sounds counter-intuitive that a 
hadron experiences more collisions in a less absorptive medium. Formally it follows from 
3), but can be explained qualitatively. For instance, if one calculated the mean number of 
collisions in a photoabsorption reaction on a nucleus using the Glauber formula, the result 
will be very small, proportional to a em . However, N co u is defined for events when inelastic 
collision took place. In this case it comes from hadronic fluctuations of the photon and is 
much larger than number of collisions given by the Glauber formulas Eq. (3). This example 
explains why N coU increases due to inelastic shadowing. 

1.2.2 How inelastic is the inelastic cross section? 

As far as one needs to calculate the deuteron-nucleus inelastic cross section it should be 
done in correspondence with the class of events selected by the trigger. The cross section 
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calculated via Glauber Monte-Carlo generator in all three experiments corresponds to the 
Glauber formula derived in Appendix A, Eq. (A. 14), where af n N is the total inelastic NN 
cross section. Then, according to derivation, Eq. (A. 14) describes the total inelastic cross 
section on a nucleus minus the part related to quasielastic nuclear excitations (with no 
hadron produced). This is not what was actually measured in any of the three experiments 
[1-3]. These experiments have different event selections and the calculations should comply 
with that. 

The STAR experiment triggers on forward neutrons from the gold [2] and detects all 
inelastic d — Au collisions including quasielastic excitation of the gold 1 . In this case, accord- 
ing to the Glauber formalism presented in Appendix A, one should rely on Eq. (A. 10) with 
a tot' = 51 m b, rather than a^ n N = 42 mb. At the same time, the two other spectrometers, 
PHENIX and PHOBOS, seem to be insensitive to large rapidity gap events, i.e. diffractive 
excitations of the deuteron and gold [1,3] 2 . Then Eq. (A. 14) should be applied with a re- 
placement af n N a^ n N — 2a^ d N — a dd N ss 30 mb, i.e. the single and double diffraction must 
be subtracted (see details in Sect. 2) [12,13]. Apparently, it makes difference whether one 
performs calculations with input cross section 51 mb, or 42 mb, or 30 mb. 

The numerator in (4), af n N is even more sensitive than the denominator to assumptions 
which inelastic channels should be included. However this does not seem to be a problem, 
since the cross section of high-p T production in pp collisions, da pp /d 2 pT, was directly mea- 
sured in the all three experiments 3 , and we consider only the nuclear modification factor 
Eq. (5) in what follows. 

1.3 The outline 

This paper is organized as follows. We present a brief and simple derivation of basic formulae 
of the Glauber model [16] in Appendix A. In Sect. 2 we treat the deuteron as a nucleus 
and generalize the Glauber model for this case. We derive formulae for the cross sections 
of different channels, perform numerical calculations and present the results in Table 1. We 
corrected the input inelastic NN cross section for diffraction and found a smaller af A than 
in [1,3], but larger than in [2]. 

Events with a tagged spectator nucleon may serve as a sensitive probe for nuclear trans- 
parency, since the spectator must propagate through the nucleus with no interaction. We 
calculate the total cross section for this channel and the transverse momentum distribu- 
tion of the spectators. On the contrary to naive expectation that noninteracting nucleons 
retain their primordial Fermi momentum distribution, we found an amazingly strong focus- 
ing effect. Namely, the nucleus acts like a lens focusing the spectators into a narrow cone 

*I appreciate the very informative communication with Carl Gagliardi on this issue. 

2 A part of diffraction might have been included into the trigger efficiency of the PHENIX spectrometer. 
Namely, double diffraction, i.e. excitation of nucleons in both the deuteron and gold can reach and fire 
sometimes the closest of the two BBC triggers covering pseudorapidity intervals r\ — ±(3 — 3.9). However, 
the main part of diffraction, single diffractive excitation of either the deuteron, or gold, hardly could reach 
the opposite hemisphere and fire both the BBC triggers simultaneously, what is the trigger condition. 1 am 
thankful to Barbara Jacak and Sasha Milov for informative and clarifying discussions of this issue. 

3 Although it is stated in [1] that the cross section da pp /d 2 pT was normalized to 42 mb, in fact it was 
measured [14, 15]. 
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with momentum transfer range of the order of the inverse nuclear radius. The transverse 
momentum spectrum of the spectators acquires typical diffraction structure having minima 
and maxima. 

Inelastic shadowing corrections are introduced in Sect. 4. First, we use the traditional 
presentation in terms of inelastic diffractive excitations in intermediate state of hadron- 
nucleus elastic amplitude (Sect. 4.1). This approach is quite restricted, being unable to 
deal with higher order scattering terms which are especially important at high energies. 
Therefore, we switch to the eigenstate representation introduced in general terms in Sect. 4.2. 
Its realization in QCD is the light-cone color-dipole approach presented in Sect. 5. 

The part of the inelastic corrections related to the lowest hadronic Fock component 
consisted only of valence quarks corresponds to diffractive excitation of resonances in usual 
terms. This contribution is analyzed and estimated numerically in Sect. 5.1. We demon- 
strate that this corrections make heavy nuclei much more transparent: instead of exponen- 
tial attenuation we found a linear dependence on the inverse nuclear thickness (Sect. 5.1.1). 
Correspondingly, we derived formulae for cross sections of different channels corrected for 
inelastic shadowing for hadron-nucleus (Sect. 5.1.2), and deuteron-nucleus (Sect. dA) colli- 
sions. In Sect. 5.3 we study possibility to improve our calculations. We tested sensitivity of 
our results to the form of the nucleon wave function, and derived formulae for the case of a 
realistic saturated dipole-nucleon cross section. 

Gluonic excitations corresponding to Fock states containing extra gluons are considered 
in Sect. 6. They correspond to diffractive excitations of large mass which are known to 
have quite a small cross section. This smallness leads to a prediction of rather weak gluonic 
shadowing, ~ 20%, and small contribution to the inelastic corrections. At the same time 
other models predict much stronger gluon shadowing (Sect. 6.2) which may substantially 
change the normalization of the dAu data. 

Since nuclear matter become more transparent due to inelastic shadowing, the number 
of participants changes as well. In Sect. 6.3 we found this effect to be sizeable. 

In Sect. 7 we sum up the effects considered so far to see how much they affect the dAu 
data. The results are presented in Table 1. We also corrected the PHENIX data for high-p T 
pions to see how important are these corrections compared to the current error bars. We 
found a considerable change: the Cronin effect for high-p^ pions disappeared. 

Our observations are summarized in Sect. 8. The main conclusion is that the current 
data for high-pr hadron production in deuteron-gold collisions are not decisive, and should 
be complemented with direct measurements of the inelastic dAu cross section. 

2 Extending the Glauber model to deuteron-nucleus 
collisions 

The basic formulae of Glauber model for hadron-nucleus collisions are presented in Ap- 
pendix A. If to treat the deuteron as a hadron, one can calculate the dA total, total elastic 
and inelastic cross section provided that the elastic dN amplitude is known. The latter can 
be calculated employing the Glauber model too. This is done in Appendix A. 
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One can do calculations differently, treating the deuteron as a system of two nucleons in- 
teracting with the nucleus. In this case one can consider more reaction channels as deuteron 
excitation etc., which are missed in the previous approach. 



2.1 The total cross section 



. We generalize Eq. (A. 5) from Appendix A for a deuteron beam as following, 




x 



1 _ f[ [i _ TP N (b-f T /2 - s k )] [l - T nN (b + f T /2 - s k )] 



fc=i 



+ 




(7) 



where tt is the transverse nucleon separation in the deuteron; |\I/d(7Y)| 2 is the deuteron 
light-cone wave function squared and integrated over relative sharing by the nucleons of 
the deuteron longitudinal momentum. It is presented in Appendix B. The effective nuclear 
thickness function, (6), convoluted with the NN elastic amplitude is introduced in (A. 6). 
We did calculations with nuclear density in the Woods-Saxon form 



with Ra = 6.38 fm and a = 0.54 fm, same as in [1] for easier comparison. The result for the 
total cross section af^ u is shown in Table 1. 

Eq. (7) is easy to interpret. The two first term in the exponent correspond to independent 
interaction of two nucleons separated by transverse distance tt- Of course, the smaller tt 
is, the stronger nucleons shadow each other, and this is accounted for by the third term. 

One can see the difference between this expression and Eq. (A. 5) (for h = d). In the 
latter cases the averaging over tt is put up into the exponent, while in the former case, 
Eq. (7), the whole exponential is averaged. We will see in Sect. 4.2 that this difference is 
a part of the Gribov's inelastic corrections, so (7) makes the first step beyond the Glauber 
approximation. 

Note that the last term in the exponent in (7) is quite small. Besides smallness of 
a ei N I ' a u>t i the exponential factor is rather small. The mean value of the exponent is 
{rj)/AB NN ~ 5. This term reduces the total d — Au cross section by 1.3% only. 

2.2 The cross section of elastic dA scattering and deuteron breakup 



Pa{t) 



3A 1 



(8) 



47ri$(l + 7r 2 a 2 /i$) 1 + exp (^) 



dA — ► pnA 
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Table 1: Results for different cross sections and numbers of collisions calculated 
using Glauber approximation (Sect. 2), corrected for inelastic shadowing related 
to valence quark fluctuations (Sect. 5.1), and for gluon shadowing (Sect. 6). The 
results including the ultimate renormalization factor K depend on the experi- 
mental set up and are different for the STAR and PHENIX experiments 





Observable 


Glauber 
model 


Valence quark 
fluctuations 


Plus gluonic 
excitations 


Correction 
factor 




<7***[mb] 


4110.1 


3701.0 


3466.2 




s 


<rf„ Au [mb] 


2422.7 


2226.6(2335.8) 


2118.3(2228.3) 




T 


Factor K in (5)-(6) 


K GL =1.04 




K Gr = 0.87(0.92) 


K=0.91(0.96) 


A 


N™ u {min.b.) 


6.9 


7.5 


7.9 




R 


afn U (tagg)[mb] 


458.4 


544.9(511.5) 


551.8(520.1) 








2.9 


4.4 


5.0 




P 


^ tt -*//[nib] 


2146.0 


1998.3(2100.1) 


1930.3(2033.7) 




H 


Factor K 


K G i = 0.92 




K Gr = 0.9(0.95) 


K=0.83(0.87) 


E 


N coU~ dlff ( min - b -) 


5.5 


5.9 


6.1 




N 




324.3 


480.2(451.5) 


498.4(470.6) 




I 


K °r diff (tagg) 


2.3 


2.9 


3.2 




X 













According to (A. 9) in order to find elastic dA cross section one should square the partial 
elastic amplitude and integrate over b, 



dA 



= d 2 b 



d 2 r T \Vd(r T ) 



1 — exp 

2 



' 2 a tot 



+ a e r^(fe)exp 



(9) 



4B NN / 

This is the square of the average of the elastic amplitude over deuteron configurations. 
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If, however, to take average of the amplitude squared, that will include also dissociation 
d — > pn, i.e. 



a d /(dA -> dA) + a d f ss (dA -> npA) = / ^ / d 2 r r |tf d (r T 

x 1 — exp 



1 
~2 



-\o%? ( TZ(U \r T ) + T»(b- \r T ) 



+ a e r^(6)exp 



2.3 The total inelastic cross section 



T 



4B 



NN . 



(10) 



Subtracting from the total cross section the elastic part one gets the cross section of all 
inelastic channels in d — A. We, however, prefer to subtract the deuteron quasielastic 
breakup too, since it is not detected by any of the RHIC experiments. Then we have, 



dA 



a 



~ - ° d d L(dA - npA) = Jd 2 bJ d 2 r T \* d (r T ) 



x 



exp 



+ 2 or 7^(6) exp 



NN 
'tot 



T^{b + \r T ) + n{b- l -r T ) 

4 ' 



AB 



NN , 



;n) 



The result of numerical calculation for this cross section is exposed in Table 1. This cross 
section covers diffractive excitations as well, therefore we place the result in the upper panel 
of the table which is supposed to be related to experiments sensitive to diffraction (STAR). 

The impact parameter distribution of the inelastic cross section is plotted by dashed 
curve in Fig. 2, and the integrated cross section is shown in Table 1. 

2.4 The cross section of nondiffractive channels 



In experiments insensitive to large rapidity gap event one should employ the inelastic 
cross section with all diffractive contribution is removed. Namely, one should also subtract 
the cross sections of quasielastic excitation of the nucleus, A — > A*, and diffractive excitation 
of colliding nucleons. 

The cross section of single (dA — > dA*) and double (dA — > pdA*) quasielastic and 
quasidiffractive nuclear excitation reads [compare with (A. 13)], 



a d J(dA - dA*) + a d q f d (dA - pnA*) = J d 2 r T \V d (r T ) 



x 



- 



i_rj[i_ r ^(fe-f T /2-4 

k=l 

i n [i-r pN (b-f T /2-s k ) 



l-T nN (b + f T /2-s k 
l-T nN (b + f T /2-s k 
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Figure 2: The impact parameter distribution of inelastic deuteron-gold 
collisions (three upper curves) including diffractive excitations (STAR 

— * 

trigger). Impact parameter b corresponds to center of gravity of the 
deuteron. The dashed curve corresponds to the Glauber approximation 
Eq. (10). The thin solid curve include inelastic shadowing related to 
excitation of the valence quark skeleton, Eq. (45). The thick solid curve 
is Gnal, it includes gluon shadowing as well. The bottom solid thick curve 
shows the difference between the Glauber and final curves. The dotted 
curve shows the range of model uncertainty and corresponds to gluon 
shadowing with R G = 03 (see Sect. 6.2). All curves are calculated with 
total cross section a™ = 51 mb. 
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AB 



NN 



n (b + + T%(b - fr) ) + <t»»T» (b) exp [ - 



T%(b + -r T )+T%(b-^r T 



NNrpN/ 



AB 



NN 



7 (r) = 1 - 
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is a correction factor hardly different from one. In what follows we do not keep the small 
terms in (13). Note that the form of Eq. (12) is analogous to that of Eq. (A. 13). 

For experiments insensitive to diffraction the quasielastic cross section Eq. (12) should 
be subtracted (11) and the result will be similar to Eq. (A. 14). However, we still miss 
the contribution of diffractive channels related to diffractive excitations of nucleons in the 
deuteron and nucleus. It is impossible to introduce consistently diffraction in the framework 
of the Glauber model which is a single channel approximation. Diffraction naturally emerges 
in the multiple coupled channel approach or in the eigenstate method introduced below. 
Meanwhile, one can use the following prescription. 

Let us expand the exponentials in (12) in small expansion parameter a^ N T^b) up to 
the first order, 
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(14) 



In order to include the possibility of diffractive excitation of nucleons in the colliding 
nuclei, one should replace in (14) 



cr. 



MN 
el 



~NN 



cr. 



el 



NN , c>~-NN NN 
a el + 2(T .<d + Cr r , 



' dd 



(15) 



in all orders of Ta(6). This is a substantial correction since at RHIC energy cry = 9 mb, 
anda e r + 2crf/ + ^/ = 21mb. 

The final Glauber model expression for the non-diffractive inelastic dA cross section 
reads, 



dA 

non—diff 



d 2 r T \V d (r T ) 



exp 



~NN 



T^b+ l -T T )+T^b- l -T T ) 



+ 4a 



~NN rpNN 
el 2 A 



(b) exp 



AB 



NN . 



(16) 



where 

~NN NN ~NN 

a in = a tot ~ °el I 17 ) 

We calculated the non-diffractive part, Eq. (16), of the inelastic d — Au cross section, 
and the result is shown in Table 1. The corresponding number of collisions also presented in 
the Table is rather small compared to the one quoted in [1]. This is mainly due to a smaller 
inelastic cross section a^L N we use. 



3 Quantum mechanics at work: illuminating focusing 
effect for spectators 

Assume that only the proton in the deuteron interacts inelastically with the nucleus, while 
the neutron is a spectator (of course all following results are symmetric relative to inter- 
change p «-> n). This is a very interesting process of simultaneous interaction and no inter- 
action. It provides direct information about nuclear transparency. Apparently, this process 
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pushes the neutron to the ultra-periphery of the nucleus where its survival probability is 
high, while the proton prefers to hit the dense area of the nucleus and interact. 

Naively, the survived spectator neutrons should maintain their primordial transverse 
momentum distribution controlled by the deuteron size. This is assumed in the Glauber 
Monte Carlo. However, quantum mechanics is at work, and the nucleus acts like a lens 
focusing spectator neutrons. The survival probability modifies the shape of the wave packet 
of the spectators in the impact parameter plane. Correspondingly, their ^-distribution 
changes. This is how elastic scattering on an absorptive target happens: it is not due 
to transparency of the target, but is caused by absorption. In the limit of a completely 
transparent target, the incoming plane wave is not disturbed and no scattering occurs. 
Absorption makes a hole in the plane wave, and one can think about the outside area of the 
incoming wave which undergoes elastic scattering. On the other hand, one can subtract the 
incoming plane wave whose Fourier transform is just a delta function (zero angle scattering) 
and the rest is a wave packet with a transverse area of the target size. A Fourier transform 
of this wave packet gives the elastic amplitude [compare with (A. 9)]. 

Thus, the spectator neutrons experience elastic scattering on the target, rather than 
simply propagate with the undisturbed primordial transverse Fermi momentum. Below we 
derive formulas which show how elastic scattering of the spectator neutrons happens, and 
perform numerical evaluation of the effect. 

We start with the cross section of this process which can be written as, 
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(18) 



where b is the impact parameter of the proton. The first factor here would be the elastic 
neutron-nucleus cross section, if it were not weighted by the second term which is the 
inelastic proton-nucleus cross section, i.e. the difference between the total and elastic and 
quasielastic cross sections (see Appendix A). 

After integration over the coordinates of bound nucleons we get, 
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(19) 
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Here we made a correction for diffract ive channels replacing a, 
valid only for those experiment which are not sensitive to diffraction (PHENIX, PHOBOS). 
Correspondingly, the numerical result for a l ^_ di jj(dA — > nX) is placed at the bottom 
part of the Table. The results at the upper part of the Table use the total elastic and 
total cross sections instead of a^ N and af n N as an input for calculations. Indeed, since 
the STAR experiment is sensitive to quasielastic nuclear excitation as well, it should be 
included, and one has to rely on Eq. (A. 13). We also demonstrate the impact parameter 
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Figure 3: The impact parameter distribution of interacting protons in 
tagged deuteron-gold collisions with spectator neutrons. The calcula- 
tion includes diffractive excitations (STAR trigger). Impact parameter 

— * 

b corresponds to the proton. The dashed curve represents the Glauber 
approximation Eq. (10). The thin solid curve include inelastic shadowing 
related to excitation of the valence quark skeleton, Eq. (45). The thick 
solid curve includes gluon shadowing as well. All curves are calculated 
with total cross section aFJ* = 51 mb. 



dependence of a°f 9 (dA — * nX) in Fig. 3. Interesting that the interacting protons in tagged 
dA collisions strongly pick at the very edge of the nucleus in spite of the large radius of the 
deuteron. This is not a trivial observation which can be probably interpreted as follows. 
The spectator neutron must be mostly outside of the nucleus. Then, for a protons which 
are close to the edge of the nucleus the interval of azimuthal angle between the proton and 
neutron is larger than for a proton deep inside the nuclear area. This phase space factor 
enhances the contribution of peripheral protons. 

To see how the spectator neutrons are distributed one can use the same Eq. (19) with 
the replacement b =>- b + r T . The result of calculations is depicted in Fig. 4. This plot 
demonstrates that the spectators have a more peripheral impact parameter distribution 
than the interacting protons, but they are amazingly close. 

We also calculated the number of collisions of the proton which underwent interaction 
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Figure 4: Tie impact parameter distribution of spectator neutrons in 
tagged deuteron-gold collisions with interacting protons. The calcula- 
tion includes diffractive excitations (STAR trigger). Impact parameter 

— * 

b corresponds to the proton. The dashed curve represents the Glauber 
approximation Eq. (10). The thin solid curve include inelastic shadowing 
related to excitation of the valence quark skeleton, Eq. (45). The thick 
solid curve includes gluon shadowing as well. All curves are calculated 
with total cross section aFJ* = 51 mb. 




in events with a tagged spectator neutron, 
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(20) 



The results for N^ff for events which include diffraction or do not are shown at the upper 
and bottom panels of Table 1 respectively. The mean value of N co u for tagged events turns 
out to be nearly a half of the minimal bias value, Eq. (3), which is for two nucleons in 
the deuteron. This contradicts the intuitive expectation that tagged events are much more 
peripheral than minimum bias inelastic collisions and the proton should have a much smaller 
number of collisions. 

To get the transverse momentum distribution of spectator neutrons, one should Fourier 
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transform the elastic neutron amplitude before squaring it, 
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(21) 



Important is that the proton inelastic interaction is incoherent, therefore we should first sum 
up coherently all amplitudes of neutron elastic scattering for the fixed impact parameter of 
the proton, then Fourier transform it, square and after all integrate over the proton impact 
parameter. This is explicitly done in (21). Apparently, integration over qr in (21) leads to 
the expression in Eq. (19). The S and D wave functions are presented in Appendix B. 
In Fig. 5 we compare the normalized differential cross section Eq. (21), 
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with the undisturbed primordial distribution of the neutron in the incoming deuteron, also 
normalized to one, 
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(23) 



The surprising observation is that the spectator neutrons have a much narrower q T - 
distribution than the Fermi motion in the deuteron. This is opposite to the usual qr- 
broadening (Cronin effect) for particles propagating though a matter [61]. In the present 
case the nucleus acts like a lens focusing neutrons. Fig. 5 also exposes quite a different 
shape of the (^-distribution of spectators having diffraction-like minima and maxima. 

If to compare the mean values of q\ of spectator neutrons with the primordial value in 
the deuteron, the difference is tremendous, about factor 20. 

(q T )s P ect = 0.00038 GeV 2 (24) 
(q 2 T ) 'deuteron = 0.0065 GeV 2 (25) 

This focusing effect is a beautiful manifestation of quantum mechanics. The intuitive 
interpretation is rather straightforward. The condition that the neutron in the deuteron 
remains intact, while the proton must interact, means that the neutron tries to pass the 
nucleus through the diluted periphery, while the proton prefers the collision to be central. 
These conflicting conditions cause a strong suppression of small-size deuteron fluctuations, 
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Figure 5: Transverse momentum distribution of spectator neutrons in 
the tagged reaction d + Au — > n + X (solid curvej, and in tie projec- 
tile deuteron (dashed curve). The inelastic reaction p + Au — > X is 
assumed to include diffraction (STAR experiment) . The calculations are 
performed in the Glauber approximation, Eq. (21). 



while large separations in the deuteron are enhanced. Apparently, such large size configu- 
rations are related to a smaller Fermi momentum and this simple observation explains the 
focusing effect. 

This explanation offers a possibility to study the correlation of the focusing effect with 
centrality of collision 4 . Suppressing 6-integration in (21) one can trace the 6-dependence 
of the focusing effect. In Fig. 6 the same comparison of two (^-distributions is shown for 
central collision b = (impact parameter of the interacting proton). The interpretation 
of central collisions is especially clear. Once the proton hits the center of the nucleus, the 
spectator neutron must be located along a rage ring outside of the nucleus with a radius 
larger than the nucleus. Correspondingly, the (^-distribution has the typical diffractive 
shape and a small width Aq T < 1/Ra- 

On the contrary to our expectations, the distributions are quite similar. The mean value 
of (<7y) = 0.00032 GeV 2 is close to our result Eq. (24) for the minimal bias sample. 

One may wonder why the minima on the (^-distribution Fig. 6 are deeper than for min- 

4 I am thankful to Alexei Denisov for this suggestion. 
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Figure 6: Same as in Fig. 5, but for a central (b = 0) proton-gold collision, 
accompanied by a spectator neutron. 



imal bias sample Fig. 5. In fact, for central collisions the minima go down to zero, since we 
neglected the real part of the elastic amplitude and the Fourier transform oscillates changing 
sign. However, the position of the minima (slightly) depend on the impact parameter of the 
collision. Therefore, when one sums up (^-distributions with different minimum positions, 
the resulting distribution will have minima which are partially filled up. 

One should be cautious comparing these predictions with data which might be contam- 
inated by non-spectator neutrons. First, the neutron calorimeters used at RHIC have a 
rather large acceptance which covers transverse momenta up to ~ 300 MeV. Therefore, 
most of the neutrons which experienced quasielastic scattering contribute as well (except 
STAR). Besides, the range of longitudinal momenta is rather large and events with diffrac- 
tive excitation on nucleons in the gold should contribute too. All such neutrons are not 
spectators and have much wider (^-distribution. 

Second, selecting central collisions in accordance with higher multiplicity, one should 
remember that central collisions are suppressed (see Fig. 3) and one should not mix them 
up with the fluctuations of multiplicity. 
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4 Inelastic shadowing corrections 



4.1 Intermediate state diffractive excitations 

The Glauber model is a single-channel approximation, it misses the possibility of diffractive 
excitation of the projectile in intermediate state illustrated in Fig. 7. These corrections 




Figure 7: Diagonal and off-diagonal diffractive multiple interactions of 
the projectile hadron in intermediate state. 



called inelastic shadowing were introduced by Gribov back in f969 [9]. The formula for the 
inelastic corrections to the total hadron-nucleus cross section was suggested in [17], 
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where 



XN with longitudinal 



(27) 



is the cross section of single diffractive dissociation hN 
momentum transfer 

_ M 2 -m\ 

This correction makes nuclei more transparent [18]. One can also see from Fig. 1 that 
(26) does a good job describing data at low energies [10, 19], since takes care of the onset of 
inelastic shadowing via phase shifts controlled by qi- Higher order off-diagonal transitions 
are neglected. Diagonal transitions (or absorption of the excited state) are important, but 
unknown. Indeed, the intermediate state X has definite mass M, but no definite size, or 
cross section. It is ad hoc fixed in (26) at cr^. It has been a long standing problem how to 
deal simultaneously with phase shifts which are controlled by the mass, and with the cross 
section which depends on the size. This problem was eventually solved in [20,21] within the 
light-cone Green function approach (see Sect. 6). 

The situation changes at the high energies of RHIC and LHC, all multiple interactions 
become important, but phase shifts vanish, substantially simplifying calculations. No exper- 
imental information, however, is available for off-diagonal diffractive amplitudes for excited 
state transitions Xi — > X 2 . A solution proposed in [22] is presented in the next Sect. 4.2. 
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There is, however, one exclusion which is free of these problems, hadron-deuteron colli- 
sions. In this case no interaction in intermediate state is possible and knowledge of diffractive 
cross section NN — * NX is sufficient to calculate the inelastic correction with no further 
assumptions. In this case Eq. (26) takes the simple form [9,23], analogous to (A. 17), 

r r dn hN 

*o2 = -2fw>frt w ± s F#). (28) 

We calculate this correction for pd collisions following [24] at s/s = 200 GeV using the 
slope B$ N = 10 GeV" 2 reduced compared to B% N = 14 GeV" 2 by 4 GeV" 2 which is the 
proton vertex contribution to the elastic slope. The upper cut off imposed by the deuteron 
formfactor on integration over M 2 is quite high at this energy and we can use the free 
diffraction cross section a^ d N = 4mb [12]. Then we find Aa^j. = — 1.75 mb 



4.2 Eigenstate method 

If a hadron were an eigenstate of interaction, i.e. could undergo only elastic scattering 
(as a shadow of inelastic channels) and no diffractive excitation was possible, the Glauber 
formula would be exact and no inelastic shadowing corrections would be needed. This simple 
observation gives a hint that one should switch from the basis of physical hadronic states to 
a new one consisted of a complete set of mutually orthogonal states which are eigenstates of 
the scattering amplitude operator. This was the driving idea of description of diffraction in 
terms of elastic amplitudes [25,26], and becomes a powerful tool for calculation of inelastic 
shadowing corrections in all orders of multiple interactions [22]. Hadronic states (including 
leptons, photons) can be decomposed into a complete set of such eigenstates \k), 

W=E*£l*>> (29) 
k 

where are hadronic wave functions in the form of Fock state decomposition. They obey 
the orthogonality conditions, 

£ K)' ^ = * hh , ; 
k 

E(*f) t$ J = (30) 
h 

We denote by f^ N = i ofj^ /2 the eigenvalues of the elastic amplitude operator / neglect- 
ing its real part. We assume that the amplitude is integrated over impact parameter, i.e. 
that the forward elastic amplitude is normalized as \f^\ 2 = Att da^ N /dt\ t=Q . We can then 
express the elastic f e i{hh) and off diagonal diffractive f s d{hh!) amplitudes as, 



/r = 2i£|*fc| o% = 2i(a)- t (31) 
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f^ N {h^h')=^Y J {0^ h k^t ■ (32) 
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Note that if all the eigen amplitudes were equal, the diffractive amplitude (32) would vanish 
due to the orthogonality relation, (30). The physical reason is obvious. If all the f^ N are 
equal, the interaction does not affect the coherence between the different eigen components 
\k) of the projectile hadron \h). Therefore, off diagonal transitions are possible only due to 
differences between the eigen amplitudes. 

If one sums up all final states in the diffractive cross section one can use the completeness 
condition (30). Then, excluding the elastic channels one gets [22,27,28], 
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As far as the main problem of Glauber approximation is the need to include off-diagonal 
transitions, one should switch to an eigenstate basis. Then each of the eigen states can ex- 
perience only elastic diffractive scatterings and the Glauber eikonal approximation becomes 
exact. Thus, all expressions for cross sections of different channels derived in Glauber ap- 
proximation in Appendix A are exact for any of the eigenstates. Then, the corresponding 
cross sections for hadron-nucleus collisions are obtained via a proper averaging of those in 
Appendix A [22,28], 
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It is interesting that the last expression for cr^ 4 is already free of diffraction contribution. 
Although only elastic and quasi-elastic cross sections were subtracted from <7^ f in Glauber 
model in Appendix A, after averaging over eigenstates it turns out that diffraction is sub- 
tracted as well. Indeed, direct averaging of the elastic cross section Eq. (A. 9) is different 
from (35) and includes coherent diffraction, hA — > XA, which cross section reads [22,28], 
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(37) 



Averaging of the quasi-elastic cross section Eq. (A. 13) leads to inclusion of diffractive exci- 
tation of the hadron h — ■> X besides excitation of the nucleus, A — > Y . 

Thus, Eq. (36) resulting from a direct averaging of the single channel inelastic cross 
section Eq. (A. 14) corresponds to the part of the total hA cross section which does not 
contain, elastic scattering, hA — * hA, coherent diffraction, hA — > XA, quasi-elastic, hA — > 
hY, and double diffraction, hA — > XY. This part of the cross section is what is measured 
as the inelastic cross section in heavy ion and p(d)A collisions at SPS and RHIC, and what 
we are going to calculate below. 

One may wonder, what is the difference between the cross sections Eqs. (34)- (36) and 
those in Glauber approximation, Eqs. (A. 5), (A. 9) and (A. 14)? The difference is obvious, 
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in the former set of equations the exponentials are averaged, while in the Glauber approxi- 
mation contains exponentials of averaged values. For instance, the total cross section in the 
Glauber approximation reads, 
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where {a tot) — r tot ■ ^ to subtract this from Eq. (34), the rest is the Gribov's inelastic 
correction calculated in all orders. Indeed, we can compare it with the expression Eq. (34) 
expanding the exponentials in (34) and (38) in multiplicity of interactions up to the lowest 
order. Employing (33) we find, 
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This result is identical to Eq. (26), if to neglect there the phase shift vanishing at high 
energies, and also to expand the exponential. 

Note that since the inelastic nuclear cross section in the form Eq. (A. 14) is correct for 
eigenstates, one may think that averaging this expression would give the correct answer. 
However, such a procedure includes possibility of excitation of the projectile and disinte- 
gration of the nucleus to nucleons, but misses possibility of diffractive excitation of bound 
nucleons which is not a small correction. We introduce a corresponding correction in the 
next section. 



5 Light-cone dipoles and inelastic shadowing 
5.1 Excitation of the valence quark skeleton 

The light-cone dipole representation in QCD was introduced in [28] where it was realized that 
color dipoles are the eigenstates of interaction and can be an effective tool for calculation of 
diffraction and nuclear shadowing. It was concluded that the key quantity of the approach, 
the cross section of the dipole- nucleon, o^ry), is a universal and flavor independent function 
which depends only on transverse separation tt and energy. Of course the energy must be 
sufficiently high to freeze variations of the dipole size during interaction, otherwise one 
should rely on the Green function approach [20,21,29] (see Sect. 6). 

This representation suggests an effective way to sum up all multi-step inelastic correc- 
tions in all orders [28]. Since dipoles are eigenstates of interaction in QCD, they are not 
subject to any diffractive excitation, and the eikonal approximation becomes exact. There- 
fore, if energy is high enough to keep the transverse size of a dipole "frozen" by Lorentz time 
dilation during propagation through the nucleus, one can write the cross sections in the form 
Eqs. (34)- (36). The averaging in this case means summing-up different Fock components 
of the hadron consisted of different numbers of quarks and gluons, and for each of them 
integration over tt (intrinsic separations), weighted with the square of the hadron light-cone 
wave function |^ / / l (rr)| 2 . We assume that the hadron does not have a "molecular" struc- 
ture, i.e. is not like a deuteron consisting of two colorless clusters. Therefore all following 
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expressions apply only to elementary hadrons. To simplify calculations, in what follows we 
rely on the quark-diquark model of the proton, neglecting the diquark size. The total cross 
section is basically insensitive to the diquark size, besides, there are many evidences that 
this size is indeed small [30,31]. 

5.1.1 Nuclear transparency 

According to the Glauber model hadrons attenuate exponentially in nuclear matter, 



where Tr, called nuclear transparency, is the survival probability of a hadron propagating 
through a nuclear matter of thickness T4. However, we know that the hadron fluctuates 
and can be viewed as a combination of Fock states of different content and size. Some of 
them having a small transverse size can easily penetrate the medium and do not attenuate 
as fast as in (40). 

Assuming that the hadronic wave function has a Gaussian form and the dipole cross 
section g(tt) oc r\ (this small-TY behavior does a good job describing hierarchy of hadronic 
cross sections and their sizes [32]) we can perform averaging in (34) and arrive at a rather 
simple expression [28], 



This explicitly demonstrates how Gribov's corrections make nuclei more transparent. Since 
exponential attenuation is much stronger than a power, for large Ta (central collisions with 
heavy nuclei) the difference might be tremendous. 

5.1.2 Cross sections 

The total cross section. 

The total hadron-nucleus cross section is modified according to (41) as, 



Although the Gribov's corrections (color transparency) make nuclei much more trans- 
parent, the modified total cross section Eq. (42) is not much smaller than the result of 
Glauber approximation Eq. (A. 5). This is because the central area of a heavy nucleus is 
"black", i.e. fully absorptive, in both cases, and the cross section is mainly related to the 
geometry of the nucleus. In other words, the exponential term in (A. 5) is very small for 
central collisions, and the total cross section is rather insensitive to even dramatic variations 
of its magnitude. 

The elastic cross section 
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The partial elastic cross section is given by the square of the averaged value of the elastic 
amplitude. We get, 
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Correspondingly, the differential elastic cross section reads, 
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The total inelastic cross section 

The cross section of all inelastic channels is given by the difference, 
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This cross section approaches the unitarity limit for T^(b) ^> 1 at the nuclear center, 
but is proportional to T\{b) at the nuclear periphery. 
Diffractive excitation of the hadron. 

The combined cross section of elastic scattering and diffraction when the hadron may 
be either excited or not, but the nucleus remains intact, is given by the average of the dA 
elastic partial amplitude squared, 
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Here we firstly averaged over the quark coordinates in the nucleons, secondly, squared the 
result, and thirdly, subtracted the elastic dA cross section [compare with (37)]. 
Diffractive excitation of the nucleus. 

The cross section of the reaction where the nucleus is diffractively excited, and the 
hadron either remains intact or is excited too, reads, 
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XY)) , (48) 



and (jgj is a cross section of single diffractive excitation of either the beam or the target; 
the double diffractive cross section a 1 ^ corresponds to diffractive excitation of both. 

Deriving Eq. (47) we made use of smallness of the elastic cross section and expanded the 
exponential. Higher orders of a^ N are neglected, but the corrections are easy to calculate. 
We also neglected the small variation of the elastic slope of the dipole-nucleon cross section 
with r T . 

Eq. (47), as one can see from (48), takes into account possibility of diffractive excitation 
of the projectile. This is a direct consequence of the eigenstate approach. In addition, we 
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also included the possibility of diffractive excitation of bound nucleons in the target. Those 
excitations are not shadowed by multiple interactions in the nucleus, since all extra particles 
produced this way stay in the nuclear fragmentation region and do not break down the large 
rapidity gap structure of the event. Therefore, they may be incorporated into a^ t N adding 
the two last terms. At yHs = 200 GeV single and double diffraction cross section are about 
equal, af d N « a^ N « 4 mb [12, 13], a% N « 9 mb, so a% N « 21 mb. 

Diffractive reactions Eq. (47)-(48) do not produce any particles at central rapidities. 
Therefore, if one wants to calculate the part of the total hadron-nucleus cross section de- 
tected experimentally, one should subtract these diffractive contributions, 

= ^-*s d A +e i(hA^XA)-a^(hA^XA*) 



qel 1 

- / ( fh "I" 1 ■'■■""> ) I _ " J et 1 l°tot l , u \ \ 

/ " U 1 , -hNmhrL\ I n , -hN rrhti.xn ( ■ \* V J 



1 + <r$ T%(b) \ + 



Since pp cross section is used as a baseline for comparison, the same subtraction should be 
done in this case too, 

(50) 
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what comes to about = 30 mb at y/s = 200 GeV. 

Then, the number of collisions at given impact parameter corrected for inelastic shad- 
owing reads, 
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5.2 Deuteron-nucleus collisions 

So far we considered the case of a colorless hadrons, but colored constituents. The specifics 
of a deuteron is that it contains two colorless clusters, nucleons. Therefore, one of the 
inelastic corrections which we already took into account in (7) is related to fluctuations 
of the deuteron size. The next step is to average over the fluctuations of the sizes of the 
nucleons. 

The total deuteron-nucleus cross section. 

Now we should average af A over the inter-nucleon separation, as well as over the nucleon 
sizes, r*i and r 2 , 

af A = 2jd%J d\ T 

where 

1 ~ [i + \^FtW+ \r T )][i + \o^n{b-\f T )\ 
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The result of calculation exposed in Table 1 is smaller than the Glauber model value. The 
difference comes from inelastic shadowing related to diffractive excitations of the colorless 
clusters in the deuteron, each consisted of three valence quarks. 

Elastic and diffractive scattering of deuterons. 
Correspondingly, the total cross section of elastic scattering and diffractive excitation of the 
deuteron has the form, 



where 
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Inelastic deuteron-nucleus collisions. 

If to subtract the elastic and diffractive cross section Eq. (54) from (52) the rest will be 
the inelastic cross section which covers all diffractive excitations of the nucleus, but not 
gold. This is what is measured at the STAR experiment. To comply with the condition of 
experiments insensitive to diffraction one should also subtract the cross section of diffractive 
excitation of the nucleus. The results read [compare with (49)], 



~dA 



a 



dA 
tot 



a 



dA 



a 



(dA 



XA) - a d J(dA 



dY) - o%(dA -> XY) 



where 



d 2 b d 2 r T ^ d (r T ) (h dA (b,r T ) 



ri,r 2 



(56) 



h dA (b,r T ) 



1 - 



[i + o$f mb + \r T )\ [i + o$? n^ -i r T )\ 

2arnjb + \r T ) 

[1 + nib+ \r T )f[l + a? ot N n (b - \t t )] 



2amib-^f T ) 

[i + n0+ \ r T )\ [i + nib - \ 



mb 



r T ) 



2a™nib)e W 



[1 + *?ot N n(b+ \r T )Y[l + ag?mb-\r T )\ 



(57) 



27 



The results of calculations of both inelastic cross sections with and without nuclear 
diffraction, as well as the corresponding numbers collisions which are rather small compared 
to what was calculated in [1], are presented in Table 1. As expected, the cross sections are 
smaller than predicted by the Glauber model, while the numbers of collisions are larger. We 
also plotted 6-dependence of af^ in Fig. 2 (thin solid curve). Comparing with the Glauber 
curve we see that this class of inelastic shadowing corrections leave the mid of nucleus 
"black" , but makes it rather transparent on the periphery. 

Production of spectator nucleons. 
Similarly, one derives an equation for the cross section of a channel with tagged spectator 
nucleons corrected for inelastic shadowing, 



*d(r T ) 
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Events with tagged nucleons are especially sensitive to the transparency of the nucleus. 
We calculated the cross section, Eq. (58), and the results, as well as the corresponding 
numbers of collisions are shown in Table 1. The effect of inelastic corrections on the impact 
parameter distribution of interacting protons in tagged events with a spectator neutron 
is demonstrated in Fig. 3. Calculation was done for inelastic proton interaction including 
diffractive excitations (STAR). As one could anticipate, the nucleus becomes much more 
transparent in the center. Indeed, for a nearly black nucleus inelastic corrections keep it 
black since transparency, or the exponential term is so small that even if it is modified by 
a large factor, the final change is very small. However, tagged events is a direct measure 
of transparency, and the inelastic corrections are maximal in this case. It is not surprising 
that N co u is quite large (considering that only one nucleon interacts. 

5.3 Towards realistic calculations 
5.3.1 Three valence quarks 

For the sake of simplicity we used so far the approximation of a quark-diquark structure 
of the proton, and neglected the diquark size. As long as the diquark is indeed as small 
as 0.2 — 0.3 fm [30,31], this approximation is rather precise even for heavy nuclei which 
hardly can resolve such a small size. However, the mean size of the isoscalar diquark is still 
a debatable issue, besides, an isovector diquark is probably a big object. Then, one may 
expect nuclear matter to be more opaque for a high energy nucleon compared to what was 
found above. 

We evaluate nuclear transparency for another extreme, i.e. for the case of a proton wave 
function symmetric in all quark coordinates, with a mean size of any diquark of the order 
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of 0.7 fm, 



|*jv(fi,f 2 ,f 3 )| 2 = - — j^ ex P — " 1 " 5(ri + f 2 + f 3 ). (59) 

\ Kr N) \ r N / 

To perform the averaging of the eikonal exponentials in (34)- (36) we need to know the 
three-body dipole cross section, which we express via the conventional qq one as, 

<73q(ri, f 2 , f 3 ) = - [(Tqqin) + <J qq {r 2 ) + <7gg(r 3 )] . (60) 

It satisfies the limiting conditions, namely, turns into cr qq (r) if one of three separations 
is zero. Assuming that cr qq (r) = Cr 2 , this cross section averaged with the wave function 
squared Eq. (59) gives a™ = Cr%/2. 

Now we can calculate the nuclear transparency averaging the eikonal exponential, 

f 3 

exp [-a 3q (ri)T A (b)}) = \{d 2 r i \^ N (f u f 2 , f 3 )| 2 exp [-<7 3(? (fi, f 2 , f 3 ) T A (b)\ 

(61) 



l + k<t N T A (b) 

We see that nuclear transparency in this case is quadratic, rather than linear function of 
the inverse nuclear thickness. For small a^f T A (b) 1 it coincides with the result of the 
quark-diquark model, Eq. (41), however falls steeper at large T A {b). This is not surprising: 
in order to make use of color transparency the whole proton has to fluctuate into a small 
transverse area, and it is more probable for a two- than three-body system. 

One can consider these results as a lower [Eq. (61)] and an upper [Eq. (41)] bound for 
nuclear transparency. We calculated different cross sections using the average of the eikonal 
exponential in the form Eq. (61) instead of Eq. (41), and the results are shown in Table 1 
in parenthesis. Unfortunately, we still do not know the proton wave function sufficiently 
well to fix this uncertainty for nuclear transparency. Nevertheless, the difference is not large 
for real nuclei. For instance, the inelastic non- diffract ive d — Au cross section presented in 
Table 1 increases by about 6%. 



5.3.2 Realistic dipole cross section 



The dipole cross section a^ q oc r\ used above is justified only for small tt, while it is 
expected to level off at large qq separations. More reliable calculations can be done using 
a realistic phenomenological cross section. A quite popular parametrization was proposed 
in [33] and fitted to HERA data for F 2 (x, Q 2 ). However, it should not be used for our 
purpose, since is unable to provide the correct energy dependence of hadronic cross sections. 
Namely, the pion-proton cross section cannot exceed 23 mb 5 . 

5 According to [34] this dipole cross section reproduced well the energy dependence of the photoabsorption 
cross section a^ t (s). This happens only due to the singularity in the light-cone wave function of the photon 
at small Tt- This is a specific property of the transverse photon wave function and is not applicable to 
hadrons. 
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A parametrization more appropriate for soft hadronic physics was proposed in [11]: 



oo(s) 



1 — exp 



T 



R 2 o(s), 



(62) 



where Rq(s) = 0.88 fm (so/s) ' 14 and s = 1000 GeV 2 . In contrast to [33] all values depend 
on energy (as it is supposed to be for soft interactions) rather than on x and the energy 
dependent parameter ao(s) is defined as, 



vo(s) = aZ(s) 1 + 



3r 2 (s) 
8 (r 2 ch ), 



(63) 



Here {r 2 ^^ = 0.44±0.01 fm 2 [35] is the mean square of the pion charge radius. Cross section 
(62) averaged with the pion wave function squared automatically reproduces the pion-proton 
cross section. The pp total cross section is also well reproduced using the quark-diquark 
approximation for the proton wave function. The parameters are adjusted to HERA data 
for the proton structure function. Agreement is quite good up to at least Q 2 ~ 10 GeV 2 
sufficient for our purposes. 

With such a dipole cross section one can perform analytic calculations expanding the 
Glauber exponentials in (34)- (37). Then the total cross section gets the form, 
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Correspondingly, the sum of elastic and diffractive deuteron scattering on the nucleus 
reads, 
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The cross section of quasielastic excitation of the nucleus with simultaneous possibility 
to excite the deuteron is given by, 
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In all these equations 
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Now one can calculate subtracting (65) and (66) from (64). However, in this pa- 
per we restrict ourselves by calculations performed above and leave this more complicated 
computation for further study. 
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6 Gluon shadowing and the triple- Pomeron diffraction 



First of all, to avoid confusion it should be emphasized that we are not talking about gluon 
shadowing in high-p-^ hadron production at xp = in d — Au collisions. This process 
exploits Bjorken x > 0.01 which are too large for gluon shadowing [4, 11]. On the contrary, 
we consider gluon shadowing in the soft inelastic d — Au collisions which are the main 
contributor to the total cross section. This process is related to much smaller x ~ 10~ 5 . 

Gluon shadowing is an important source of inelastic corrections at very high energies. It 
is pretty clear if one employs Eq. (26). The part of the diffraction which corresponds to the 
triple-Regge graph IPIPM, or the lowest order Fock component consisted only of valence 
quarks, has a steep M-dependence, da^ jdM 2 oc 1/M 3 . Therefore the integral over M 2 in 
(26) well converges, the minimal momentum transfer vanishes at high energies, and this 
part of inelastic corrections saturates. 

The triple-Pomeron (IP IP IP) part of diffraction which corresponds to the Fock state 
containing at least one gluon, is divergent at large masses, da^ jdM 2 oc 1/M 2 , since the 
gluon is a vector particle. The cut off is imposed by the nuclear formfactor in (26), i.e. the 
condition q^ < 1/Ra- As a result of the divergence, this part of the inelastic corrections 
rises as ln(s/s ) and reaches a substantial value at the energy of RHIC. 

Eikonalization of the lowest Fock state \3q) of the proton done in (34)- (36) corresponds to 
the Bethe-Heitler regime of gluon radiation. Indeed, gluon bremsstrahlung is responsible for 
the rising energy dependence of the phenomenological cross section (62), and in the eikonal 
form (34)- (36) one assumes that the whole spectrum of gluons is multiply radiated. However, 
the Landau-Pomeranchuk-Migdal (LPM) effect [36, 37] is known to suppress radiation in 
multiple interactions. Since the main part of the inelastic cross section at high energies is 
related to gluon radiation, the LPM effect becomes a suppression of the cross section. This 
is a quantum-mechanical interference phenomenon and it is a part of the suppression called 
Gribov's inelastic shadowing. The way how it is taken into account in the QCD dipole 
picture is inclusion of higher Fock states, \3qG), etc. Each of these dipoles is of course 
colorless and its elastic amplitude on a nucleon is subject to eikonalization. 

As already mentioned, Eq. (26) should not be used at high energies as it misses all higher 
order multiple off-diagonal transitions, and incorrectly (ad hoc) calculates diagonal ones. On 
the other hand, the eigenstate expressions Eqs. (34)- (36) is not safe to use either. Indeed, 
the significant part of the integral over M 2 in (26), next to the upper cut off, corresponds to 
a finite qi- In other words, the fluctuation valence quarks + gluons is not frozen by Lorentz 
time dilation during propagation through the nucleus. 

6.1 The Green function for glue-glue dipoles 

A proper treatment of a quark-gluon fluctuation "breathing" during propagation through 
a nucleus is offered by the light-cone Green function formalism. In this approach the ab- 
sorption cross section, as well as the phase shifts, are functions of longitudinal coordinate. 
This is also a parameter-free description, all the unknowns are fixed by the comparison with 
other data. We employ this approach and calculate gluon shadowing following Ref. [11]. 
The key point which affects the further calculations is the nonperturbative light-cone 
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wave function of the quark-gluon Fock state, 




(68) 



Here we assume (as usual) that the gluon is carrying a negligible fraction, a G < 1, of 
the quark momentum. This wave function is quite different from the perturbative one 
which is the same as in light-cone description of the Drell-Yan process [11,38]. The latter, 
employed for calculation of diffractive gluon radiation (the triple- Pomeron term), results in 
overestimation of data for large mass diffraction more than by order of magnitude. This 
problem has been known since the 70s as the puzzle of smallness of the triple-Pomeron 
coupling. The way out is to make a natural assumption that the parent light-front quark 
and gluon experience a nonperturbative interaction which squeezes that quark-gluon wave 
packet, and therefore reduces the dipole cross section. The parameter r in (68) controls 
the strength of the real part of the light-cone potential which is chosen in a Gaussian form. 
Fit to diffractive data pp — > pX leads to the value of the mean transverse q — G separation 



y (r 2 ) = r = 0.3 fm. This conclusions goes along with the results of nonperturbative 
models, like the instanton vacuum model [39], and lattice calculations [40], which found a 
similar small size for gluonic fluctuations. Such a semihard scale, 1/ro also leads to quite a 
steep energy behavior of the radiation cross section and well explains data for the total and 
differential elastic cross sections of pp scattering [41]. 

Apparently, smallness of r leads to quite a weak shadowing for Fock states containing 
gluons. As a consequence, we expect rather weak gluon shadowing, what is not a surprise 
in view of the close connection between diffraction and shadowing. As long as the gluon 
clouds around valence quarks small, the Gribov's corrections are suppressed. Besides the 
fluctuations containing gluons become heavy and the onset of gluon saturation takes place 
at much smaller x than usually expected. 

The mean quark-gluon separation r ~ 0.3 fm is much smaller than the quark separation 
in light hadrons. For this reason one can neglect the interferences between the amplitudes of 
gluon radiation by different valence quarks. Since the gluon contribution to the cross section 
corresponds to the difference between the amplitudes of \qqqG) and \qqq) components, the 
spectator quarks cancel out. Then the radiation cross section is controlled by the quark- 
gluon wave function and color octet (GG) dipole cross section. 

Thus, the contribution to the total hadron-nucleus cross section which comes from gluon 
radiation has the form, 
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where «g is the fraction of the quark momentum carried by the gluon; 
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Here the energy and Bjorken x are related as s = 2m^E q = 4/xTq. 

The second term in (70) corresponds to the triple-Pomeron part of the inelastic correction 
Eq. (26) written in impact parameter representation. The amplitude of diffractive gluon 
radiation qN — > GqN is proportional to ty qG (r, a G )a GG (r). A glue-glue dipole emerges in 
this expression because this is not elastic scattering, but a production process. Its amplitude 
comes from the difference of the scattering amplitudes of different Fock components of the 
quark [38], \q) and \qG), which is a dipole cross section of a color octet-octet dipole, qq — G. 
Since the size of the qq pair is irrelevant for gluon shadowing, we neglect it and replace 
the qq by a gluon (see in [11,42]). Therefore the second term in (70) can be interpreted as 
production of a qG pair at the point z±, then propagation of this pair with varying transverse 
separation up to point z 2 where it converts back to the quark. 

Propagation of the dipole of varying size through the absorptive medium between points 
z\ and z 2 is described by the Green function G GG (f 2 , z 2 ; ri, z±). It satisfies the two dimen- 
sional Schrodinger equation, 



i~T— G GG (f 2 ,z 2 ;f 1 ,z 1 ) 
dz 2 



A(f2) +V(f 2 ,z 2 ) 



G GG {r 2 ,z 2 \r x ,z x ) , (71) 



2E q a G (l-a G ) 

where imaginary part of the light-cone potential is related to absorption in the medium, 

Im V(f, z) = ~ a GG (r) p A (b, z) . (72) 

For further calculations we assume that the quark energy is E q = s/6mjv, but the results 
are hardly sensitive to this approximation. 

Perturbative calculations treating a quark-gluon fluctuation as free particles overesti- 
mates the cross section of diffractive gluon radiation (or the triple-Pomeron coupling) more 
than by an order of magnitude. The only way to suppress this cross section is to reduce the 
mean transverse size of the fluctuation. This is done in [11] via introduction of a real part 
of the light-cone potential in (71), 

2 

Re V(r, z) = T -^— l , (73) 

1 J 2E q r%a G (l-a G ) 1 ; 

where parameter tq = was fitted to data for single diffraction pp — > pX. 

The gluonic dipole cross section a GG (r,s) is assume to be different from the qq one, 
Eq. (62), only by the Cazimir factor 9/4. To simplify calculations we rely on the small-r 
approximation, a GG (r, s) ~ C GG (s) r 2 , where C GG (s) = dcr GG (r, s)/dr 2 =Q . This approxima- 
tion for the dipole cross section is justified by the small value of Tq « 0.1 fm 2 . 

In the case of a constant nuclear density, Pa(j') = Pa ®{Ra — r ), the solution of Eq. (71) 
has a form, 
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where 
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Integrations in (70) can be performed analytically, 

P(a G , b) = Re \n(W) , 
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The first term in (70) is a part of the nuclear cross section calculated in Bethe-Heitler 
limit, i.e. without gluonic inelastic shadowing. Therefore it is included in the nuclear cross 
sections calculated so far. The new inelastic shadowing correction comes from the second 
term in (70). Its fraction of the total pA cross section is depicted in Fig. 8. The onset of 
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Figure 8: Ratio of the gluonic inelastic shadowing correction (minimal 
bias) to the total nuclear cross section as function of cm. energy y/s. 



shadowing is delayed up to yfs ~ 20 GeV. We believe that this result is trustable since the 
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Green function approach treats phase shifts and attenuation in nuclear matter consistently. 
Nevertheless, in order to get an idea about the scale of theoretical uncertainty we also 
evaluated the magnitude of gluon shadowing using the known values of the triple-Pomeron 
coupling and equation (26). The results are quite similar, in both cases the gluon shadowing 
correction is pretty small [11], ~ 20% at the energy of RHIC. Such a weak shadowing is 
a direct result of smallness of the parameter tq = 0.3 fm which we use. This seems to be 
the only way to suppress diffractive gluon radiation corresponding to the triple-Pomeron 
contribution, and to reach agreement with data on diffractive dissociation pp — > pX. For 
this reason, all effects related to gluons, including saturation, or color-glass condensate, are 
quite suppressed. 

Naturally, the inelastic correction in (70), (77) varies with impact parameter vanishing 
on the very periphery and reaching a maximum at central collisions. At small 7a (6) the 
inelastic correction is proportional to T\(b) while the partial amplitude is proportional 
to TaQ)). Therefore, the ratio linearly rises with T A {b) (see in [43,44]) with a coefficient 
approximately equal to 0.2 fm 2 . For very large T\(b) the correction may even exceed the 
rest of the cross section, then apparently higher order corrections must be added to stop 
this growth. Such a saturation is not important for real nuclei, therefore we use the linear 
parametrization Ro(b) = 1 — A^(6) = 1 — 0.2T A (b) for further calculations. 

The valence quark part of the inelastic shadowing corrections makes the nucleus more 
transparent, i.e. reduces the elastic scattering amplitude as one can explicitly see comparing 
the corrected amplitude Eq. (42) with the Glauber form Eq. (A. 5). However, both approach 
the black disc limit for large T A (b)a^ ^> 1. Important question is whether this is still true 
after inclusion of gluonic corrections. 

Eq. (70) has the typical form of a nonlinear equation like Glibov-Levin-Ryskin evolutions 
equation (GLR) [45], or in the dipole form Balitsky-Kovchegov equation (BK) [46]. The 
second term on the right-hand side of (70) corresponds to glue-glue fusion in GRL or the 
multiple interaction in the nucleus in BK equations. We calculated the correction in the 
lowest order using the uncorrected dipole cross section o" GG (r), i.e. the undisturbed free 
gluon density. Next iterations would be to implement the corrected gluon density (at larger 
x, however), or accij'), into the second term in (70). This procedure leads to the BK 
equation which solution is still a challenge. However, due to smallness of the correction, 
20%, we do not expect large higher order corrections and the saturated solution should not 
be very different from our result which we employ in further applications. 

On the other hand, if gluon shadowing emerging from the first order iteration is very 
strong, like it was found in [47-49], it should be substantially reduced by next iterations 
which effectively play role of self-screening. Namely, as long as the gluon density is reduced 
at small x, one cannot use in Eq. (26) the cross section of diffractive dissociation on a free 
nucleon target. It is suppressed by the same gluon shadowing (at larger x though). The 
stronger is the gluon shadowing the more important is this self-screening effect. It was 
missed in calculations [47-49] which grossly over-predicted the strength of gluon shadowing. 

Now we are in a position to correct our previous calculations for the gluonic part of 
inelastic shadowing which we fix at 20%. We do it replacing ^(r^) =>- Rc{b) &qq{ r T)) 
where Rcip) = 1 — A in (b) is the suppression factor related to gluon shadowing. This simple 
prescription is based on the intuitive expectation that a dipole interacts with a less number of 
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gluons in the nucleus than the eikonal model assumes. Indeed, for small separations ty, the 
dipole cross section reads [50], (Jg q (r T ) = (iT 2 /3)a s r 2 ^G(x,r T ), i.e. it is indeed proportional 
to the gluon density which is reduced in nuclei. More motivations for this procedure can be 
found in [43,44]. 

The results for nuclear cross sections corrected for gluon shadowing are shown in Table 1 
and depicted in Fig. 2-4 by thick solid curves. 

6.2 More models for gluon shadowing 

Although we predict quite a modest qluon shadowing effect and therefore a rather small 
inelastic shadowing correction, many models predict much stronger effects. One can call it 
theoretical uncertainty, if one treats all models equally (though some of them are probably 
more equal than others [51]). 

For instance, the popular event generator HIJING contains a Q 2 -independent gluon 
shadowing [52] which is a factor 0.3 at x ~ 10 -5 . With such a dramatic gluon shadowing 
we get the impact parameter dependence of the inelastic cross section depicted by dotted 
curve in Fig. 2. The corresponding correction factor K = 0.65 for the PHENIX data. 

If to treat shadowing in terms of the dipole approach, it is clear that shadowing is a 
monotonic function of Q 2 , since the size of the dipole can only rise towards the soft limit. 
This is confirmed by DGLAP evolution of nuclear shadowing in the perturbative domain. 
Therefore, one can use gluon shadowing predicted by different models at the starting scale 
Q of the order of 1 — 2 GeV 2 as a bottom bound for the shadowing correction expected in 
the soft limit. We calculate Bjorken x for the RHIC energy y/s = 200 GeV and Q 2 = 1 GeV 2 . 

A strong gluon shadowing was predicted in [47,48], R G = 0.3 — 0.4. HERA data for 
diffraction 7*p — > Xp was used as an input in equation (26) modified for 7*A collisions. 
The statistics of this data is much lower than in proton diffraction pp — > pX and not 
sufficient for reliable determination of the triple-Pomeron coupling. Different solutions for 
this coupling fitted to DIS diffractive data vary dramatically [53]. Besides, as is mentioned 
above, the gluon self-screening missed in [47,48] should significantly reduce the effect of 
gluon shadowing. 

Explicit calculations of gluon shadowing via gluon dipoles was performed in [49]. The 
gluon shadowing corresponding to the RHIC energies was found at Rg ~ 0.6 what leads to 
correction factor K = 0.78 for the Phenix ratio Rha- This calculation, however, also does 
not include the gluon self screening, and is based on the assumption that gluon and quark 
dipoles have identical distribution functions. 

A strong gluon suppression was also found in a model with an early onset of strong 
saturation [54] which characteristic scale is a steep function of energy, Q 2 oc (l/rr) - 252 . 
It was assumed in the KLM approach [6] that for Q 2 < Q 2 S gluon density xGa(x, Q 2 ) 
is proportional to Q 2 R\ with a factor which was taken from the McLerran-Venugopalan 
model [55] at x — 10 _1 . Such an oversimplified picture exhibits a strong gluon shadowing. 
If to compare the xGa{x, Q 2 ) with the GRV parametrization [56] at x ~ 10 -4 it turns out 
to be strongly suppressed by factor Rq = 0.42. In this case the correction factor in (5) is 
K = 0.72. 

Such a diversity of model predictions suggests a conclusion that the current data for 
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deuteron-gold collisions [1-3] cannot resolve in a model independent way the dilemma, 
whether final state interaction or initial conditions is the main source of hadron suppression 
in heavy ion collisions. Indeed, if the latter were true, it would unavoidably lead to a 
substantial reduction of af n A and the ratio Eq. (5) (compared to the Glauber model). 



6.3 Number of participants 

Although the concept of number of participants originates from a naive treatment of mul- 
tiparticle production called wounded nucleon model, it is a widely used characteristics of 
centrality of collisions. We are not going to dispute here its meaning, but just to see how it 
is affected by the inelastic corrections 6 relying on its formal definition, 



dN p (s,b) 



d 2 b 



Gl 



T A {s- b) {l - exp[-<r™ T B (b)} } + T B (s) {l - exp^a™ T A {s - b)] } 



(79) 

where s is the parameter of collision of nuclei A, B. We use subscript Gl to emphasize that 
it corresponds to this model which is inspired by the Glauber model (although they have 
nothing in common). 

Apparently, inelastic shadowing corrections should reduce N p since nuclear matter be- 
comes more transparent. The corrected expression for N p reads, 



dN p (s, b) 
cPb 



= <7„-, 



X 



T A {s-b)T B {b) 
R G (b) 

1 + R G (b) ag N T B (b) ' l + R G (s-b)a? n N T A (s-b) 



+ 



R G (s-b) 



(80) 



Here the gluon shadowing factor R G is function of impart parameter according to calcula- 
tions in [11,43,44] and the parametrization used above. We present the correction to the 
"Glauber" expression defines as, 



J shad 



(b) 



Rg(P) aZ N T B (b) 
1 + Ra(b) < N T B (b) 



exp 



—a. 



NN 



T B (b)}} , 



(81) 



in Fig. 9 depicted by solid curve. As one could expect, the correction factor peaks at the 
nuclear periphery and approaches one at large impact parameters. 

If to compare with Glauber calculations employing the incorrect inelastic cross section 
af n N = 42 mb, the correction is even larger, as is demonstrated by dashed curve in Fig. 9. 



7 Cronin effect: renormalizing the data 

Cronin effect for high-p T pions at ^Js = 200 GeV was predicted in [4] to be a rather small 
enhancement, about 10% at the maximum. The smallness of the effect is due to the change 
of the mechanism of high-p^ particle production which takes place at the RHIC energies. 

6 I am thankful to Larry McLerran suggested to look at this parameter. 



37 




Figure 9: Solid line is the correction factor Eq. (81) for inelastic shad- 
owing to the number of participants in p — Au collisions as function of 
impact parameter. Dashed curve also includes a correction to a^ n N (see 
text). 



At lower energies (SPS, CERN) different bound nucleons contribute to this hard process 
incoherently. The nuclear enhancement is due to initial/final state pr-broadening of partons 
propagating through the nucleus. This broadening should not be translated into a modi- 
fication of the parton distribution in the nucleus since /c^-factorization is broken [57]. At 
high energy an incoming light-cone fluctuation which contains a high-p^ parton is freed via 
coherent interaction with many nucleons in the target. It turns out that such a coherent 
mechanism leads to a weaker Cronin enhancement than the incoherent one. This is why 
calculations [58,59] missing this effect of coherence, predict a stronger Cronin effect. 

The PHENIX data for neutral pions [1] are depicted in Fig. 10 by full points in compar- 
ison with the predicted ratio [4]. However, as it was stressed above, the normalization of 
the data is based on Glauber model calculations which are subject to different corrections, 
all of which have negative sign. As a result, the data should be renormalized according to 
Table 1 by multiplying the experimental values by coefficient K = 0.83. The corrected data 
are shown by open circles. 

Cronin effect on a deuteron 

Theoretical predictions has been done so far for pA collisions. In order to compare 
models with dA data one should make sure that the Cronin enhancement on the deuteron 
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Figure 10: The Cronin ratio RAu/diPr) for pions. Open circles show the 
results of PHENIX with normalization base upon Glauber model cal- 
culations of the inelastic d — Au cross section using af n N = A2mb [1]. 
Full points show the same data corrected for a proper value of inelastic 
NN cross section and Gribov's inelastic shadowing. The error bars in- 
cludes statistic and systematic uncertainties. The curve is the prediction 
from [4]. 

itself is a small corrections. We evaluate the ratio, 

RM - 2d^m~ T ' (82) 

in the limit of short coherence length which gives an upper estimate for the effect. Following 
[4] the pd cross section at high p T is given by the following convolution, 

&pd(PT) = J2 fyp ® F i/d ® °ij-+kl ® D h/k , (83) 
i,j,k,l 

where Fi/ p and FjM are the distributions of parton species i,j dependent on Bjorken 
and transverse momenta of partons in the colliding proton and deuteron respectively. The 
beam parton distribution Ff is modified by the transverse momentum broadening of the 
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projectile parton due to interaction with another nucleon in the deuteron. The broadening 
of the mean transverse momentum squared reads [60,61], 



A (k 2 T ) 



da m (r T ) 



dr% 



(T) 



(84) 



r T =0 



where (T) is the mean nuclear (deuteron) thickness covered by the projectile parton before 
or after the hard collision, 



(T) 



n hN 
a tot 



d 2 s ReT hN (s) 



*d(0) 



(85) 



where we neglected the elastic slope B^n compared to the nuclear radius squared. For the 
parton distribution functions in a nucleon we use the leading order GRV parametrization 
[56]. 

We calculated this ratio using the computer code for the Cronin effect developed in [4] 7 , 
and the deuteron wave function ^(rr) described in Appendix B. 

The results for R d / p (p T ) are depicted in Fig. 11. Indeed, the Cronin enhancement is only 
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Figure 11: Cronin ratio Rpdipr) calculated at ^Js = 200 GeV using the 
formalism developed in [4]. 



2%, and can be neglected comparing d — Au data with predictions done for p — Au. 



8 Summary and conclusions 

The main observations and results of this paper are: 

7 I am thankful to Jan Nemchik who performed this calculation for pd collisions. 
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The current normalization of inclusive high-p^ cross section in deuteron-gold collisions 
measured at RHIC is based on Glauber model calculations of the inelastic dAu cross 
section which is subject to Gribov's inelastic shadowing corrections. Importance of 
these corrections is not debatable, they have solid theoretical ground and are confirmed 
by precise measurements [10,19] (see Fig. 1). These corrections, Eq. (26), have negative 
sign, i.e. make nuclear medium more transparent, and they rise with energy. 

First of all, the Glauber calculations must be improved. The inelastic iViV cross section 
used as an input should be corrected for diffraction. For experiments insensitive to 
diffraction (PHENIX, PHOBOS) the cross section should be reduced from af n N = 
42 mb down to a^ n N ~ 30 mb. On the contrary, if an experimental trigger detects 
diffraction (STAR) this cross section should be increased up to a™ = 51 nib. This 
modification results in a correction factor Kqi presented in Table 1. 

There are two types of inelastic shadowing corrections. One corresponds to diffractive 
excitation of the valence quark skeleton, or nucleonic resonances, and is related to the 
1P1P1R triple-Regge graph. We calculated this correction, Eqs. (42), (61), (64), using 
the light-cone dipole representation which effectively sums up all orders of multiple 
interactions. 

Another type of inelastic shadowing is related to diffractive gluon bremsstrahlung, 
or to the soft limit of gluon shadowing in nuclei related to the IPIP1P triple-Pomeron 
diffraction. We performed calculations in Sect. 6 using the solution for the Green func- 
tion, Eq. (74), describing propagation of a glue-glue dipole through nuclear medium 
and found a rather weak gluon shadowing for gold, about 20%. At the same time, 
other models predict much stronger gluon shadowing ranging up to corrections of 70% 
(Sect. 6.2). 

Altogether, we expect a reduction of inelastic dAu cross section compared to what 
was used for normalization of high-pr data at RHIC. We conclude that the published 
data should be corrected by factor K which is about 0.8 for PHENIX and about 0.9 
for STAR (see Table 1). The renormalized data for pions do not possess any more the 
Cronin enhancement. This correction factor might be even smaller, down to 0.65 if to 
use a stronger gluon shadowing predicted by other models. 

One should admit that current data for high-p^ hadron production in dAu collisions 
at RHIC cannot exclude in a model independent way the possibility of initial state 
suppression suggested in [6], although that would contradict the author's personal 
viewpoints. 

Probably the only way to settle this uncertainty is a direct measurement of either the 
cross section of high-pr pion production in dAu collisions, or the inelastic dAu cross 
sections at RHIC. 

A very sensitive test of models for inelastic shadowing offer tagged events with a 
spectator nucleon. In situation when direct measurement of dAu inelastic cross section 
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is difficult, this might be a way to restrict models and narrow the band of theoretical 
uncertainty. The relative fraction of these events 20% measured in [2] create apparent 
problems for models with strong gluon shadowing which predict a much larger fraction. 
Even with our weak shadowing this fraction ranges between 23% and 26%. However, 
one should make it sure that the detected neutrons are really spectators, what is not 
the case currently (see discussion in Sect. 3). 

• We found a beautiful quantum- mechanical effect: the nucleus acts like a lens focusing 
spectators. In spite of the naive anticipation that nucleons which escaped interaction 
retain their primordial Fermi momentum distribution, there is a strong narrowing 
effect substantially reducing the transverse momenta of the spectators. Besides, the 
distribution acquires the typical diffract ive maxima and minima (see Figs. 5-6). 
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9 Appendix 

A Glauber model glossary 

The hA elastic amplitude at impact parameter b has the eikonal form, 



where {sj, Zj} denote the coordinates of the target nucleon Nj. iT hN is the elastic scattering 
amplitude on a nucleon normalized as, 



A 



v hA (b; or,, Zj }) = i - n [i - r hN (b - s k ) 



(A.l) 



k=l 




(A.2) 
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A.l Heavy nuclei 



In the approximation of single particle nuclear density one can calculate a matrix element 
between the nuclear ground states. 



0=1 



oo 

-Jd 2 sY hN (s) I dzp A (b-s,z) 



where 



pa$i,z 1 ) = /n^iM{^})i 2 , 



(A.3) 



(A.4) 



is the nuclear single particle density. 

Total cross section. The result Eq. (A.3) is related via unitarity to the total hA cross 
section, 



hA 
'tot 



2Re J d 2 b |l - 
2 j d 2 b |l -exp 



l-^j d 2 sT hN {s)T A {b-s) 



-2^(1-^ PP )T» 



(A.5) 



where p pp is the ratio of the real to imaginary parts of the forward pp elastic amplitude; 

2 



T h A (b) 



n hN 
"tot 



d 2 s Rer h "(s)T A (b-s) ; 



(A.6) 



and 



/oo 
dzp A (b,z), (A.7) 
-oo 

is the nuclear thickness function. We use exponential form of T hN (s) throughout the paper, 
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exp 



2B 
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(A.8) 



where B hN is the slope of the differential hN elastic cross section. Note that the accuracy 
of the optical approximation in (A.5) is quite high for gold, ~ 1CT 3 , so we use it throughout 
the paper. We also neglect the real part of the elastic amplitude in what follows, since it 
gives a vanishing correction ~ p 2 j A 2 I Z . 

Elastic cross section. As far as the partial elastic amplitude is known, the elastic cross 
section reads, 

r 1 

1 — exp 



hA 
T el 



d'b 



(A.9) 



Total inelastic cross section. Apparently it is given by the difference between the 
total and elastic cross sections, 



®in = o% - a h el A = Jd 2 b{l- exp [-o% T h A (b)] } . 



(A.10) 
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This includes all inelastic channels when either the hadron or the nucleus (or both) are 
broken up. 

Quasielastic cross section. As a result of the collision the nucleus can be excited 
to a state |F). Summing over final states of the nucleus and applying the condition of 
completeness, one gets the quasielastic cross section, 



= d 2 b 
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Here we extracted the cross section of elastic scattering when the nucleus remains intact. 
Then in the fist term of this expression we make use of the relation, 
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and arrive at, 
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(A.12) 



(A.13) 



Inelastic nondiffractive cross section. If one is interested in the fraction of the total 
inelastic cross section (A. 10) which covers only reactions with production of new particles, 
one should exclude the nucleus break up to nucleons and nuclear fragments. That is the 
quasielastic cross section, Eq. (A.13), 



hA 
prod 



Jd 2 b{l-e W [-off T>)]} 



(A.14) 



This additional subtraction makes sense only for experiments which miss the non-production 
break up of the nucleus. If, however, all inelastic events are detected, including diffractive 
(production and non-production channels) excitations of the nucleus (check with [2]) one 
should rely on Eq. (A. 10) for the inelastic nuclear cross section. 

Diffractive cross section. One needs to know this cross section in order to subtract 
it also from the inelastic cross section, since diffractive events escape registration at p(d)A 
collisions at SPS and RHIC. The Glauber approximation is valid only for a single channel 
problem. One can extend it to include diffraction properly introducing phase shifts due to 
longitudinal momentum transfer. However, one needs to know the cross section of interaction 
of the produced diffractive excitation with nucleons. This goes beyond the reach of the 
Glauber model, and instead of further ad hoc development of the model, we solve this 
problem within the eigenstate method in Section 4.2. 



A. 2 Proton-deuteron collisions 

Apparently, Eq. (A. 5) should not be applied to light nuclei, in particular to a deuteron. 
Instead one should use, 

o& = 2o™ + Aof* , (A.15) 
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where 

Ao& = -2 J d 2 b J d 2 r T \^ d (r T )\ 2 T hN (b + f T /2) T hN (b- r T /2) . (A.16) 

One can switch via Fourier transform to momentum representation in each of these three 
factors and perform integration over fx and b. The result has a form of a one-dimensional 
integral [16], 

9 r rlrr NN 

Aag = --j d 2 q T F d (4q 2 ) , (A.17) 

where F d (q 2 ) is the charge formfactor of the deuteron. We neglected the correction ~ 1CT 3 
due to the nonzero real part of the forward AW amplitude. Note that s in (A.17) is the 
deuteron diameter, rather than the radius. This is why the formfactor argument is kq\. 
We use parametrization of the deuteron formfactor from [24], 



F d (q 2 ) = 0.55 e- aq T + 0.45 e - ^ , (A. 18) 

and a = 19.66 GeV" 2 , (3 = 4.67 GeV~ 2 . 

Using a^J* = 51 mb at y/s = 200 GeV, and the elastic slope B NN = 14 GeV -2 we found 
the total pd cross section, af^ = 97 mb with Glauber correction Agj = —5 mb. Since at this 
point a correct proton-deuteron cross section is needed, we have to go beyond the Glauber 
approximation and add the inelastic correction considered in Sect. 4. We show that it is 
equivalent to adding the differential cross section of single diffraction, pN — > XN, to the 
elastic one in (A.17). This increases the value of the shadowing correction by 1.75 mb. 
Finally, we arrive at the cross sections, 

o** = 95.15 mb 

at = a? t--^- = 84.9mb. (A.19) 

167T B pd 

Interesting, the inelastic cross section is not affected by the Glauber correction, it is even 
slightly larger that the sum of two inelastic AW cross sections. The slope from the differ- 
ential elastic pd cross section was measured and fitted in [62] , 

a<J el = \ a tot) e B pd t+C pd t 2 / A 

dt 16tt ' y ' ' 

where 

B pd = b + b ± In s pd , (A.21) 

with parameters b = 32.8 ± 0.6 ( GeV 2 ) and &i = 1.01 ± 0.09 ( GeV 2 ). Parameter C pd = 

54.0 ± 0.9 ( GeV~2) was found to be energy independent. At the energy of RHIC B pd = 

44.1 GeV~ 2 , and we use this value in (A.19). 

The inelastic cross section Eq. (A.19) contains inelastic diffractive channels like quasielas- 
tic break up of the deuteron, pd — > ppn, and excitation of the nucleons pd — > Xd, pd — > pY, 
and pd — > XY. For the experiments insensitive to diffraction (PHENIX,PHOBOS) those 
channels must be subtracted. 
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B Deuteron wave function at rest and Lorentz boosted 



To perform calculations for interaction of a high energy deuteron, one should not use the 
three-dimensional deuteron wave function, but needs to know the light-cone deuteron wave 
function expressed in Lorentz invariant variables, the transverse n — p separation tt and 
the light-cone fraction a = Pn/pt °f the deuteron momentum carried by a nucleon. One 
cannot get this wave function by a simple Lorentz boost from the rest frame of the deuteron, 
where the 3-dimensional wave function is supposed to be known, to the infinite momentum 
frame. Deuteron is not a classical system, under a Lorentz boost it acquire new constituents 
which are quantum fluctuations. These constituents build up higher Fock components. This 
makes the procedure of Lorentz boost extremely complicated. There is, however, a practical 
recipe suggested in [63] and widely accepted. To the best knowledge of the author, it works 
rather well for nonrelativistic systems (nuclei [64], heavy quarkonia [65], etc.) 

The idea is straightforward, first, to express the deuteron wave function in momentum 
representation, 

M<1) = J dPre^MV , (B.l) 

via the light-cone variables in the rest frame of the deuteron. To do it one should connect 
the three-dimensional momentum squared with the effective mass of the cc pair, q 2 = 
M 2 /4 — m 2 N , expressed in terms of light-cone variables 

In order to change integration variable to the light-cone a one uses their relation, q^ = 
(a — l/2)M(q T ,a), and get a Jacobian which can be attributed to the definition of the 
light-cone wave function, 

(2 , 2 \3/4 
{.Qt "t m N) 

Applying this procedure to the S and D-wave radial wave functions one gets, 

^ =► U(f T ,a); 
r 

^1 W(f T ,a) . (B.4) 
r 

This dependence on a is important for exclusive final states, for instance deuteron disso- 
ciation to nucleons with definite longitudinal momenta. However, for most of applications 
in this paper we need to know the r T - retribution integrated over a, 

i 

\^d(r T )\ 2 = J da[u 2 (r T ,a) + W 2 (r T ,a)} , (B.5) 
o 

The result of this is identical to the simple integration over longitudinal variable in the rest 
frame of the nucleus, 



— oo 
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We use the contemporary deuteron wave functions which employ the Nijmegen-93 po- 
tential [66] 8 . 
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